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We investigate the properties of an effective Hamiltonian with competing interactions involving 
spin and chirality variables, relevant for the description of the trimerized version of the spin-1/2 
kagome antiferromagnet. Using classical Monte Carlo simulations, we show that remarkable behav- 
iors develop at very low temperatures. Through an order by disorder mechanism, the low-energy 
states are characterized by a dynamical freezing of the chiralities, which decouples the lattice into 
"dimers" and "triangles" of antiferromagnetically coupled spins. Under the presence of an exter- 
nal magnetic field, the particular topology of the chiralities induces a very slow spin dynamics, 
reminiscent of what happens in ordinary spin glasses. 

PACS numbers: 75.10.Nr, 75.10.Jm, 75.10.Hk 



I. INTRODUCTION 

The understanding and the characterization of the 
physical properties of frustrated magnetic systems is a 
central issue in modern theory of condensed matter. In- 
deed, in the last years, it has become evident that the 
magnetic interactions can give rise to a large variety of 
different unconventional phenomena, and that they can 
also be related to other phases of matter, like, for in- 
stance, superconductivity. 

From an experimental point of view, an im- 
pressive number of materials with unconventional 
magnetic properties has been discovered, including 
quasi-one-dimcnsional and ladder systems^*^ quasi-two- 
dimensional compounds with highly frustrated magnetic 
interactions j 3 - systems of weakly coupled dimers in var- 
ious geometries and fully three-dimensional frus- 
trated systems based on corner sharing tetrahedra£ 

A particular role in the list of frustrated mag- 
netic systems is played by the magnctoplumbitc 
SrCr8_ x Ga4+ x Oi9, in which the anomalous magnetic 
properties are commonly attributed to planes of Cr 3+ 
ions on the so-called kagome lattice, a lattice which 
can be described as a triangular lattice of triangles (see 
Fig. In this compound, although a small quantity 
of disorder is present in the Chromium planes, there are 
evidences that the physical properties are not much af- 
fected by the presence of impurities. Nonetheless, linear- 
susceptibility measurements did not find any sign of con- 
ventional magnetic order down to T « AK^ despite a sig- 
nificant antiferromagnetic exchange coupling J w 100/v 
between Cr 3+ ions, deduced from the Curie temperature. 
On the other hand, important insight into a possible spin- 
glass phase at very low temperatures comes from the non- 
linear susceptibility ^3 measurements by Ramirez and 
collaborators^ In usual spin-glass systems, it is known 
that X3 1 ~ [T ~ ^g) 7 : where T g is the spin-glass freez- 
ing temperature and 7 is typically between 1 and 
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In Ref. [HI it was found that %3 diverges at T g « 3.3-ftT, 
with 7 w 2, indicating a possible freezing transition at 
this temperature. It is worth noting that the divergence 



of X3 a t a finite temperature indicates the importance 
of three-dimensional effects, since two-dimensional spin 
glasses are not expected to order at finite temperature^ 

From a theoretical point of view, the possibility of a 
spin-glass behavior in the spin-1/2 Heisenberg model on 
the kagome lattice was suggested by Chandra and col- 
laborators^ The kagome lattice exhibits both frustra- 
tion and a low coordinance number, and, therefore, the 
Heisenberg antiferromagnet on the kagome lattice is a 
good candidate to find unconventional magnetic behav- 
iors at very low temperature. The most problematic issue 
in the kagome antiferromagnet comes from the infinite 
number of classically degenerate ground states^ As a 
consequence, the linear spin-wave spectrum of magnetic 
excitations possesses a whole branch of zero modest 
Moreover, exact diagonalizations on finite-size systems 
have shown that the spin-spin correlations decrease very 
rapidly with distance>i£iiLi& In addition, the series ex- 
pansion from the Ising limit and high-temperature scries 
point to the absence of magnetic order whereas, ac- 
cording to large- N approaches, the ground state should 
be disordered with unbroken symmetry^! 

The huge degeneracy of the classical ground states for 
the kagome antiferromagnet opens the possibility of the 
so-called order by disorder effect^ Indeed, it is known 
that both thermal and quantum fluctuations can lift clas- 
sical degeneracies in favor of particular states, which 
break some symmetry. Examples of this phenomenon are 
the ordering effect in the antiferromagnetic Ising model 
on the f.c.c. lattice due to thermal fluctuations^ or the 
establishment of a collinear order in the J1—J2 model 
due to quantum and thermal fluctuations^ 3 - 

In particular, thermal fluctuations could be effective 
to partially lift the classical ground state degeneracy by 
favoring states with largest entropy. In highly frustrated 
antifcrromagnets, like the kagome lattice, the resulting 
energy landscape could still be very complicated, with 
many different "valleys" separated by high energy bar- 
riers, possibly giving the system a spin-glass behavior. 
In ordinary spin-glass materials, the disordered nature of 
the system gives rise to frustrated magnetic interactions 
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FIG. 1: The trimerized kagome lattice. The solid and dashed 
lines indicate the antiferromagnetic coupling J and J', re- 
spectively. The numbers 1, 2 and 3 indicate the site indexing 
inside the elementary triangles which defines the gauge (see 
section [HJ . 



and, therefore, to many topologically different configura- 
tions that are energetically equivalent but separated by 
large (infinite) free energy barriers. The resulting physi- 
cal properties show peculiar behaviors, with a very slow 
spin dynamics 

For the classical kagome model with full rotation sym- 
metry in spin space, no evidence in favor of a spin-glass 
behavior could be found^S& and it seems that one has to 
go away from this very minimal model to have a chance to 
get a spin-glass-like behavior. Since SrCr 8 _ x Ga 4+x Oi9 is 
made of spins 3/2, it seems natural to turn to the quan- 
tum version of the model. In this paper, we consider an 
approximate way of including quantum fluctuations into 
a classical description. The starting point is the effective 
Hamiltonian obtained in Ref. [^3 for the spin- 1/2 Heisen- 
bcrg model on a trimerized kagome lattice (see Fig. |y, 
which contains both spin and chirality variables. This is 
a modified version of the spin-1/2 Heisenberg model on 
the kagome lattice, in which the exchange couplings take 
two different values J and J' according to the pattern of 
Fig-d It is worth noting that this is actually the relevant 
case for the description of SrCr8_ x Ga4 +x Oig, since the 
presence of a triangular lattice between pairs of kagome 
layers leads to two kinds of bond lengths. 

In Ref. I27I by using an appropriate mean- field approx- 
imation, this effective Hamiltonian has been useful to 
explain a number of important properties of the low- 
lying excitations of the kagome lattice, in particular the 
increase of the low-energy singlet states with the num- 
ber of sites N like 1.15 . Unfortunately, this effective 
Hamiltonian is far from being fully investigated because 
the interactions couple the spins and the chiralities in 
a non-trivial way. Nevertheless, given the very compli- 
cated nature of the problem, even at the classical level, 
the interplay between spins and chiralities can give rise 
to an unconventional behavior. In this respect, for any 
calculation in the true quantum case, it is very important 



to know the underlying classical states, and, therefore, a 
first fundamental step is to establish the classical picture. 

Recently, an interesting step toward the comprehen- 
sion of the interplay between two different kinds of classi- 
cal variables has been achieved by using a simple Ashkin- 
Teller model on a square latticed It was pointed out 
that the presence of two local degrees of freedom might 
be a key element to prevent any kind of long-range or- 
der and to obtain a very slow spin dynamics (reminis- 
cent of a spin-glass-likc behavior) at low-enough tem- 
perature. Hopefully, this finding can be seen as a first 
and crude explanation of the low-temperature behavior 
of SrCrg_ x Ga4+ x Oi9, which does not show any magneti- 
zation down to very low temperature, but presents some 
promising evidence for a spin-glass behavior. It is also 
worth noting that, although in SrCrs- x Ga4+ x Oi9 the 
Cr 3+ ions have spin S = 3/2, we consider the case treated 
in Ref.|H suitable for S = 1/2,22 

In order to go beyond the simple model considered in 
Ref. I28L based on Ising variables, and to further clarify 
the possibility to obtain a spin-glass behavior in a more 
realistic system that does not contain explicit disorder 
but only frustration between two distinct local variables, 
we investigate the effective model of Ref. [2?] in the clas- 
sical limit. 

The paper is organized as follows: In section ^ wc 
introduce the model, in section ITTT1 we present the results 
for the case of zero external magnetic field, in section llVl 
we present the results of the magnetization in the pres- 
ence of an external magnetic field, and in section Ivl we 
discuss the results and give our conclusions. 

II. THE MODEL 

Our starting point is the effective model obtained 
in Ref. |53 for the quantum spin-1/2 antiferromagnetic 
Heisenberg model on the trimerized kagome lattice. In 
this modified version of the usual spin-1/2 kagome an- 
tiferromagnet, the lattice is decomposed such as to em- 
phasize the role of the elementary triangles building the 
kagome lattice. Indeed, the spins are coupled through 
two different constants J and J' as depicted in Fig. ^ 
In the limit J' / J <C 1, the lattice can be seen as a tri- 
angular lattice of weakly coupled triangles and, by using 
ordinary perturbation theory, it is possible to derive an 
effective Hamiltonian in the subspace of the ground states 
of the triangles. The four-fold degenerate ground state 
of a triangle can be described by two spin-1/2 degrees 
of freedom, a spin a and a chirality f. With these vari- 
ables, the effective Hamiltonian is defined on a triangular 
lattice of N sites and readsi^i 

Hf = J'Y^ n Wiii (!) 
HJj = hl-2(a ij rr+alrt)}[l-2(0 ij rr+0f j rf)}, 



where denotes the sum over pairs of nearest- 

neighbor triangles. The complex parameters ay and /3y 
depend on the type of bond: They take the values 1, u) or 
tu 2 [tu = exp(27ri/3)] when the original spins in triangles 
i and j involved in the bond (i,j) sit at site 1, 2 or 3 
respectively, with the convention of Fig. ^ Note that we 
have the freedom to choose the numbering of the sites 
within each triangle, and that the convention of Fig. ^ is 
different from that of Ref. lU 

With a little algebra, we write the effective Hamilto- 
nian in a more compact way. Indeed, the operators t 4 + 
and t~ can be rewritten by introducing if and rf: 

r+ = Tf +i7f, 

Hence, a i3 T~ + a^rf = (ay + a^)rf + i(a| - ay)rf , 
and therefore: 

[2rf if ay I 

o^ T 4 + a % T t = \ -if ~ V3if if ay = to 2 
{ .;• x 3ir;' if ay =w. 

The notation is further simplified by introducing three 
unitary vectors e\ = (1,0), e*2 = (— |,— ^), e*3 = 

I 2' 2 



2t5 


■ ej 


if 


«ij 


= 1 


2t, 


■ e 2 


if 


ay- 




2t, 


' e 3 


if 


ay 


= LO. 



With this notation, the Hamiltonian reads: 

n o S = If E' ■ ^ ~ ' - 4% • 7*-), (2) 

where the vectors ey are to be chosen among e\, e*2, and 
e*3. In this Hamiltonian, we treat and as classical 
vectors with norm 1 and 1/2, respectively. By adopting 
the gauge depicted in Fig.^ we end up with a triangular 
lattice having each bond characterized by a vector e^, 
with fj, = 1,2,3. In Fig. El the different lines indicate 
which e M has to be taken for ey in the effective Hamilto- 
nian. 

Hereafter, we consider the classical limit of the effective 
Hamiltonian and we prefer to normalize both spin and 
chirality to unity, which leads to 

(id) 

In this Hamiltonian, the presence of the chiralities has the 
effect of changing the coupling between spins. We will see 
that this leads to a complicated interplay between spins 
and chiralities and gives rise to a very rich physics in the 
low-temperature regime. 
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FIG. 2: Triangular lattice on which the effective Hamiltonian 
is denned. The unitary vector for the bond is indicated by 
solid lines (e M = ei), dashed lines (e M = S2), and dotted lines 
(e M = S3). 



III. ZERO-FIELD PROPERTIES 

In this section we present the results of our model when 
no external magnetic field is applied. We provide ev- 
idences that, for low-enough temperatures, the system 
undergoes a partial order by disorder transition to a state 
where the chiralities arc essentially frozen in a given con- 
figuration. 



A. Ground state manifold 

In order to investigate the interplay between spins and 
chiralities, let us first study the low-temperature prop- 
erties of the system. This is achieved by performing 
classical Monte Carlo simulations using a heat-bath al- 
gorithmic In the context of our model, the heat-bath 
algorithm is efficient down to very low temperatures 
(T/J' « 10~ 6 ), and, therefore, it is possible to have in- 
sight into very low-energy phases. The details of this 
algorithm are discussed in Appendix^ For most calcu- 
lations, we considered an 18 x 18 triangular lattice with 
periodic boundary conditions, and we checked the con- 
clusions on larger sizes. However, as we will show in the 
following, for all temperatures, the physical properties 
are governed by short-range correlations and, therefore, 
the 18 x 18 lattice is enough to obtain insight into the 
thermodynamic limit. 

The first important outcome is that the energy per site 
converges to e = — J'/2 when T — * 0. The particularly 
simple value of the ground state energy deserves some 
further investigation. If we look back at the Hamilto- 
nian defined by Eq. we can find several classes of 
states that have energy e = —J'/2. For instance, when 
all the chiralities are parallel to a given vector 7 in the 
XY plane defined by the vectors e^, then the two factors 
involving the chiralities in Eq. (j2J) become equal, leading 
to a non-negative coupling between the spins, and the 
energy is minimized by a 120° Neel state for the spin 
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vectors. Then, the energy for an TV-site system becomes: 
E =-^ N E (l-27-^) 2 = -y/V, (4) 

^1=1,2,3 

independently of the direction 7 of the chiralities in the 
XY plane. Therefore, we obtain a large class of states 
which have, beside the usual rotational symmetry of the 
Neel state of the spins <?, an extra continuous symmetry 
in f space. 

Another important class of states is defined by those 
configurations that have chiralities anti-parallel to the 
vectors e M , independently on each site. In this case, the 
coupling of one site to four out of its six neighbors van- 
ishes. Indeed, we have that is equal to e*i for the bond 
with two neighbors, e*2 with two other neighbors and e*3 
for the last two bonds (see Fig. [2J. For any choice of 
Tfe among the three possible — e M , the effective magnetic 
coupling vanishes for the four bonds where ¥k ■ efej = ^ 
since it contains the factor (1 — 2&kj ■ t%)(1 — 2e^j • Tj). 
Moreover, having taken fx- = — e„, the effective magnetic 
coupling with the remaining two neighbors is non-zero 
only if Tj is equal to — e„ = ffc. In this case, the effective 
interaction between the two sites is equal to J'. 

A given site k can therefore interact with its neigh- 
bors in three different ways. It might be coupled to two 
neighbors, I and m, in which case these two neighbors 
will also be coupled, because they have parallel chiral- 
ities and ei m = e^i = &km- In that way, we create a 
coupled triangle, which is completely disconnected from 
the rest of the lattice. The site k might also be coupled 
to one neighbor only, forming a system of two coupled 
spins. The third possibility is to have the site completely 
disconnected from its neighbors. A typical coupling pat- 
tern is shown in Fig. |3| Hereafter, we denote by dimer 
the system of two coupled spins. It is worth mention- 
ing that because of the classical nature of the spins, the 
dimer is a classical object, with gapless excitations, and 
it has nothing to do with the quantum dimers, which are 
gapped objects. 

When the chiralities are opposite to the e^, the cou- 
pling pattern will only exhibit isolated sites, dimers and 
triangles. Moreover, any of these patterns that one can 
draw on the triangular lattice can also be achieved by 
choosing the appropriate chiralities. Having fixed the 
chiralities, the energy is minimized with opposite spins 
on the dimers and the energy per spin is e = — J'/2; on 
an isolated triangle, the sum of the spins needs to be zero, 
the energy per spin being again e = — J'/2; finally iso- 
lated sites do not contribute to the energy. Hence, among 
the possible configurations of chiralities opposite to e^, 
those that do not have isolated spins have the ground 
state energy per site e = —J'/2. 

A direct inspection of the chiralities and spins in our 
Monte Carlo simulations shows that when the system 
is cooled down abruptly from infinite temperature, that 
is by starting at the chosen temperature from a random 
configuration for the spins and the chiralities, states with 




FIG. 3: Generic coupling pattern for random chiralities op- 
posite to the vectors e M . Notice that, because of the periodic 
boundary conditions, some of the bonds (indicated by dashed 
lines) connect sites on opposite sides of the lattice. 

no obvious pattern in the chiralities and energy per site 
e = —J'/2 are achieved. Therefore, the two particu- 
lar classes of states described above are only two among 
many other ground states. On the contrary, by cooling 
down the system through a finite number of tempera- 
tures and letting it thermalize at every step before cooling 
down to the next nearby temperature, we obtain a config- 
uration with disconnected triangles and dimers. Within 
the heat-bath Monte Carlo, this fact strongly indicates 
the presence of a partial order by disorder effect: Among 
all the possible ground states, the states with large en- 
tropy are selected by thermal fluctuations. In the follow- 
ing, we will describe in more details the entropic selection 
that appears at low-enough temperature, describing how 
this transition affects the spin dynamics. 

B. Entropic selection of a family of ground states 

At temperature T/J' = 10~ 6 , the Monte Carlo results 
show that the chiralities are frozen on each site in the 
XY plane opposite to the vectors in such a way as 
to create a collection of disconnected dimers and trian- 
gles. The selection of this specific class of states can be 
explained by its large entropy. Within each decoupled 
triangle and dimer, the spins are totally free to rotate 
as long as their sum vanishes independently of the spins 
of the other triangles and dimers. Thus, all those states 
have a huge entropy at low-enough temperature. It can 
easily be proved that both for dimers and triangles, every 
site carries one degree of freedom that does not change 
the energy. This point is discussed in more details below 
in the section devoted to the specific heat. As the system 
is cooled down slowly, these states with largest entropy 
are selected. The ground state degeneracy is partially 
lifted by thermal fluctuations and an entropic order by 
disorder effect is observed. 

In Fig. 01 we show how the non-zero couplings are typ- 
ically distributed over the lattice at T/J' = 1CP 6 . We do 
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FIG. 4: Typical coupling pattern at T/J' = 1(T 6 : The solid 
lines between the sites indicate a non-zero coupling between 
the spins. Notice that, because of the periodic boundary con- 
ditions, some of the bonds (indicated by dashed lines) connect 
sites on opposite sides of the lattice. 



not have any evidence that the system prefers to create 
triangles or dimers: We have performed several simula- 
tions starting from different random orientations, and in 
all cases we have reached a state which contains both 
triangles and dimers. The average number of triangles 
is 17 on an 18 x 18 site cluster, and the distribution is 
peaked around this number with a mean square displace- 
ment of 3. We suspect that this simply corresponds to 
the distribution for all coverings of the triangular lattice 
by triangles and dimers. In any case, this observation 
agrees with the fact that, at low temperature, triangles 
and dimers have the same entropy per site. 

To have more insight into this entropic selection, we re- 
port the chiralities in the XY plane for decreasing tem- 
peratures in Fig. [S] As the system is cooled down be- 
low T* ss 0.05J', the chiralities start to display a short- 
range order in the XY plane, getting more and more 
anti-parallel to the e M vectors. As far as we can tell, this 
cross-over temperature does not depend on the cooling 
rate, provided it is not too fast. The results reported in 
Fig. [S] have been obtained after performing 10 4 sweeps 
over the lattice at each intermediate temperature start- 
ing from a random configuration, but the same behaviors 
have been observed with different numbers of sweeps. 

Since the chirality configurations that can be selected 
by the order by disorder mechanism form an infinite but 
discrete set, one can expect that, once such a state has 
been reached, the system will somehow be trapped in this 
state. This is indeed the case on short time scales. How- 
ever, whether this corresponds to a true freezing transi- 
tion requires a careful study of the low temperature dy- 
namics. For instance, in 2D models of spin glasses with 
short range interactions, the actual freezing on arbitrarily 
large time scales only takes place at zero temperature. So 
we now turn to a careful analysis of the low-temperature 
spin and chirality dynamics. 
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FIG. 5: Configuration of the chiralities in the XY plane for 
different temperatures. 
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FIG. 6: Spin-spin (upper panel) and chirality-chirality (lower 
panel) correlations at T/J' = 1CP 3 for two non-nearest- 
neighbor sites with parallel (full symbols) and non-parallel 
(open symbols) chiralities. 



C. Low-temperature spin and chirality dynamics 

Below T* , spins and chiralities behave very differently. 
Whereas the chiralities are frozen on short time scales, 
spins are expected to continue to have a fast dynam- 
ics since, once the lattice is decoupled into disconnected 
dimers and triangles, the spins belonging to each of these 
elementary building blocks are free to rotate as long as 
the total spin of the dimcr or of the triangle remains 
zero. Spin-spin and chirality-chirality correlation func- 
tions confirm this fact. In Fig. |SJ we show typical non- 
nearest- neighbor correlations <?, (t) • Bj (t) and n (t) ■ fj (t) 
as a function of the time after the system has reached 
equilibrium at T/J' = 10~ 3 . Throughout, the time unit 
is one Monte Carlo sweep over the entire lattice. In- 
deed, the spin-spin correlation function for non-nearest- 
neighbors fluctuates wildly on short time scales, while the 
chirality-chirality correlation function is essentially con- 
stant. By contrast, the spin-spin correlations for a pair 
of coupled nearest-neighbor sites are fixed to —1 or —1/2 
for dimers or triangles, respectively (not shown). 

On short time scales, the chiralities are frozen below 
T* and although there is a huge class of energetically and 
cntropically equivalent states, once one configuration is 
achieved, the system cannot switch to another one. The 
different configurations must therefore be separated by 
large energy barriers and the chiralities remain trapped 
in a small region of their phase space. A similar behavior 
is observed in ordinary spin-glass materials, where the 
spin degrees of freedom are completely frozen below a 
given temperature and different states are separated by 
infinite (in the thermodynamic limit) energy barriers. In 
this case, the transition temperature can be characterized 
by the divergence of the spin-glass susceptibility 
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1 • 



where the average over the time (Monte Carlo steps) is 
done over a sufficiently large period t e xp, corresponding 
to the time scale of the experiment. In ordinary spin 
glasses, Xsg * s zer0 ( verv small on finite systems) for 
temperatures above the critical temperature, where the 
spins have a fast dynamics, and becomes finite below the 
critical temperature, indicating that the spins are frozen 
in a given configuration. In particular, at low tempera- 
ture, Xsg g rows linearly with the number of sites N. 

By analogy, in our model, we can define a spin-glass 
susceptibility for the chirality pseudo-spins: 



Xsg 



ij i cxp *=i 



(6) 



(5) 



and try to identify a transition temperature from 
the high-temperature disordered phase to the low- 
temperature frozen phase. It turns out that, like in 2D 
disordered spin models, the freezing takes place at finite 
temperatures only for finite time scales, and that the ac- 
tual phase transition, if any, presumably takes place at 
zero temperature. To see this, we have plotted in Fig. 
the chirality susceptibility \sg as a function of tempera- 
ture for different experimental times i e xp, corresponding 
in this case, to the time we have waited at each temper- 
ature during the cooling process. For a given time scale, 
one can clearly identify two regimes: At high tempera- 
tures Xsg/-^ i s verv small, due to the dynamical char- 
acter of the chirality variables, whereas for T/J' small 
enough, there is a clear enhancement of Xsg/^' m di- 
cating a frozen dynamics of the chiralities. However, the 
temperature at which this change of behavior takes place 
depends on the time of the simulation. To quantify this 
dependence, we have plotted the freezing temperature Tf 
at which the susceptibility reaches half the value it will 
take at zero temperature as a function of i cxp . Clearly 
Tf decreases with i exp , and for times i e x P > 10 5 , the 
two are related by £ oxp = to exp(Eo/Tf), or equivalently 
T f = E /log(t exp /t ), with E ~ 0.07 J' and t ~ 6500. 

This fit shows that, at any temperature, the chirality 
configuration will change on large enough time scales, 
and that the process is thermally activated with an ac- 
tivation energy E . This energy can thus be interpreted 
as the typical free-energy barrier between two configura- 
tions selected by the order by disorder mechanism. As 
in the model of Ref[2^, the freezing temperature Tf has 
only a weak logarithmic dependence on the time scale. 

Note that these results arc not affected by the size 
(18 x 18) used in the simulation. Indeed we have checked 
in a number of cases that the conclusions are qualitatively 
and quantitatively similar for 27 x 27 and 36 x 36 clusters. 
In particular, for a given t exp , the temperature Tf does 
not change upon increasing the size of the cluster. 

To summarize, we found that spins and chiralities have 
very different low temperature dynamics in zero field. At 
a cross-over temperature T* ~ 0.05J', chirality configu- 
rations with maximum entropy are selected. Below this 
temperature, the chirality dynamics is activated with a 
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FIG. 7: Chiral spin-glass susceptibility per site Xsg/N as 
a function of the temperature for different time scales i exp . 
Inset: Variation of the freezing temperature Tj with t exp . 
Lines are guides to the eye. 
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FIG. 8: Specific heat CV per spin as a function of the temper- 
ature for different lattices: 18 x 18 (circles), 27 x 27 (squares), 
and 36 x 36 (triangles). Lines are guides to the eye. 



an energy barrier Eq ~ 0.07 J' and it is already extremely 
slow at T — 0.01J', while the spins retain a very fast dy- 
namics down to the lowest temperatures. The spins do 
not behave as free spins however, as we shall see in the 
following section, and the freezing of the chiralities has 
important consequences on the response of the spins in 
the presence of an external magnetic field. 



D. Specific heat 

On very general grounds, we can expect that the 
cross-over from the high-temperature disordered phase to 
the low-temperature phase where chirality configurations 
have been selected by the order by disorder mechanism 
is marked by a peak in the specific heat. Whereas for 
a true phase transition the specific heat diverges at the 
critical temperature, for freezing phase transitions, like 
in spin-glass systems, the specific heat exhibits a rather 
broad peak near the freezing temperature^ This is in- 
deed what we observe in our model. We calculate the spe- 
cific heat using C v = ((E 2 ) - (E) 2 )/T 2 , where E = E(T) 
is the internal energy. The specific heat per spin shows a 
broad peak which does not increase with the size of the 
lattice for T ~ T* ~ 0, 05 J' (sec Fig.®. 

The low-temperature value of the specific heat, namely 
3/2 per site, is a very nice confirmation of the picture. 
For classical variables, each quadratic mode gives a con- 
tribution to the specific heat equal to 1/2 in units of fee. 
Thus, if both spin and chirality were ordered at low tem- 
perature, we would expect the zero temperature specific 
heat per site to be equal to 2. According to our pic- 
ture however, only the chiralities are ordered. In a given 
configuration that minimizes the entropy, each chirality 
can fluctuate around its reference direction — e M with two 
quadratic modes corresponding to two orthogonal direc- 
tions, and its contribution to the specific heat is equal to 



1. However, within a triangle or a dimer, the spins are 
free to rotate as long as their sum is equal to zero. This 
means that they have zero-energy modes which do not 
contribute to the specific heat. For a dimer, the com- 
mon direction of the vectors can rotate freely, while the 
relative direction gives rise to two quadratic modes, so 
one effectively gets one quadratic mode per spin. In a 
triangle, if we consider a configuration where the sum of 
the spins is equal to zero, the orientation of a given spin 
is arbitrary, which gives rise to two zero-energy modes, 
and the others are free to rotate around the direction of 
the first one, which gives a third zero-energy mode. So 
we are left with 3 quadratic modes, i.e. one quadratic 
mode per spin, like in a dimer. So the spin contribution 
to the specific heat is 1/2 per site, and the total specific 
heat per site is equal to 3/2. Note that the entropy is 
thus the same for triangles and dimers, as stated above. 

The reduction of the spin contribution to the specific 
heat by a factor 2 is quite remarkable. This should 
be compared for instance with the classical Heisenbcrg 
model on the kagome lattice, in which case the zero tem- 
perature specific heat is reduced from 1 to 11/12 (see 
Ref|H. 



IV. FINITE-FIELD PROPERTIES 

In the absence of an external magnetic field the chiral- 
ities show a behavior reminiscent of spin glasses. In or- 
dinary spin-glass materials, the freezing of the degrees of 
freedom, together with a complicated energy landscape, 
has consequences on the magnetic properties of the sys- 
tem. Even though in our model the spins retain a fast 
dynamics even below T*, the freezing of the chiralities 
might still induce non-trivial properties when an external 
magnetic field is present. In the following, we consider 
an external magnetic field that is coupled to the spins 
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T/J' 



FIG. 9: The equilibrium magnetization per spin in units of h 
as a function of the temperature, for h/J' = 10 -2 (squares), 
2 x 10 -2 (triangles), 5 x 10~ 2 (crosses), 10 _1 (circles). Lines 
are guides to the eye. 

and we study the dynamical response of the spins. 

A. Equilibrium magnetization 

The Hamiltonian in the presence of an external mag- 
netic field h along the z axis reads: 

H eS = Hf + h'£'*l. (7) 

i 

First, we concentrate on the equilibrium magnetization 
of the system as a function of the temperature. In or- 
der to find the equilibrium magnetization per site m, we 
perform a very long simulation, cooling down the sys- 
tem through a finite number of intermediate tempera- 
tures in the presence of the external magnetic field. At 
each temperature, we equilibrate the system, by consid- 
ering a large number of Monte Carlo steps: This kind 
of protocol is known as a field-cooled (FC) experiment. 
Within this procedure, we are able to obtain accurate re- 
sults for the equilibrium curve for the magnetization as 
a function of the temperature. 

In Fig. we report the results for the magnetization 
per site as a function of the temperature, for different val- 
ues of the external magnetic field h. We can distinguish 
three different regimes. In the first one, common to all 
the cases and corresponding to T/J' > 0.1, the chiralities 
are randomly distributed and are free to rotate, retaining 
a quite fast dynamics. Therefore, (e*^ • f)) sa 0, implying 
that the effective Hamiltonian reduces to a Heisenberg 
model on a triangular lattice, with an effective antiferro- 
magnctic coupling J cS ps J'/9 [see Eq. ©]. 

In the second regime, whose temperature range 
strongly depends on the value of the magnetic field, 
the magnetization decreases and shows a plateau whose 
width depends on the actual value of the external mag- 
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FIG. 10: Configuration of the chiralities in the XY plane for 
T/J' — 10~ 6 and an external magnetic field h/J' = 1CP 2 . 

netic field. This behavior can be explained by the estab- 
lishment of a short-range order of the chiralities, which 
become opposite to the vectors e^, like in the h = case. 
In this regime, the free energy is minimized because of 
the high entropy of the spins, the physical properties are 
mostly determined by disconnected dimers and triangles 
and the actual value of the magnetization depends on the 
number of dimers and triangles, being m = h/(2J') for 
the former case and m — h/ (3 J') for the latter case. 

By further lowering the temperature, the entropic part 
of the free energy becomes less important and the chiral- 
ities rearrange by slightly modifying their ordering, al- 
lowing the spins to increase their magnetization as the 
system reaches its actual ground state for finite h. As 
T — > 0, the magnetization per site goes to m max = h/ J' . 
A direct inspection of the chiralities, see Fig. ^| shows 
that their low-temperature configuration is only slightly 
modified from the h — case. However, this is suffi- 
cient to change considerably the magnetization. The low 
temperature T m {h) at which the magnetization finally 
increases indicates that there is a very small energy dif- 
ference between the actual ground state and the states 
described by dimers and triangles. 

As shown in Fig. 13 the qualitative behavior of the equi- 
librium magnetization is the same for different values of 
h. However, as the external field gets bigger, the magne- 
tization starts increasing towards its saturation value for 
higher temperatures. Therefore, T m (h) strongly depends 
upon the external magnetic field and increases with h. 
On the contrary, the formation of dimers and triangles, 
determined by the maximum of the magnetization at 
T/J' ~ 0.1, does not depend on the external magnetic 
field, indicating that T* does not change with the ap- 
plication of h. The important fact is that, at a given 
temperature below T* , the magnetization is no longer a 
linear function of the applied external field, as typical 
magnetization curves show in Fig. ^2 

For small magnetic fields, the system is in the regime 
with disconnected dimers and triangle, which have a lin- 
ear magnetization m = hj (2 J') and m = h/ (3J 7 ), respec- 
tively. Therefore, the total magnetization is also linear in 
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FIG. 11: Magnetization per site as a function of the external 
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h with a coefficient that depends on the actual number 
of triangles and dimers. For larger h, the magnetization 
goes through a highly non-linear region, characterized by 
a complicated interplay between spins and chiralities. In 
this regime, the entropy is maximized with chiralities op- 
posite to the vectors e^, whereas the energy is minimized 
with configurations that lead to a bigger magnetization. 
Eventually, the chiralities order in such a way as to al- 
low the spins to have a magnetization m = h/J' and 
the system is again in a linear regime. A consequence of 
this competition between different states is visible in the 
dynamical properties of the system. 



B. Dynamical magnetization and slow spin 
dynamics 

In order to investigate the interplay between spins and 
chiralities and its effect on the dynamical properties of 
the system, we perform a zero-field cooled (ZFC) numer- 
ical experiment: The system is slowly cooled down to 
the temperature of interest in the absence of an exter- 
nal magnetic field. Then, an external magnetic field is 
turned on and we observe the magnetization as a func- 
tion of time. Fig. 1121 shows the result for three different 
external fields h/J' = 0.02, h/J' = 0.07 and h/J' = 0.15 
at T/J' = 10~ 3 . At equilibrium, the first field will bring 
the system into a regime with disconnected dimers and 
triangles and magnetization h/(3 J') <m< h/(2J') (cor- 
responding to the plateau region of Fig. |SJ. The third 
field will generate a state with chiralities slightly dif- 
ferent from the opposite vectors and magnetization 
m ~ h/J' (corresponding to the bump of Fig. [5] at very 
low temperatures), whereas the second field represents 
the intermediate regime. 

For h/J' = 0.02 we observe that the system acquires 
its final magnetization in a very small number of Monte 



Carlo steps, through a very fast process, typical of non- 
frustrated systems. This can be understood because, just 
before turning on the magnetic field, the system already 
displays a collection of dimers and triangles. When the 
external field is turned on, the spins quickly arrange to 
minimize the energy. In this process, the chiralities do 
not change from their previous equilibrium position. On 
the contrary, when a field h/J' = 0.07 or h/J' = 0.15 
is switched on, the final equilibrium magnetization is 
reached only after a much longer simulation time. As 
in the previous case, the system very quickly acquires a 
magnetization corresponding to the magnetization of a 
system of triangles and dimers. However, after this first 
very rapid process, a much slower dynamics takes place, 
during which the chiralities and spins find a compromise 
to minimize the free energy. Indeed, as soon as the chi- 
ralities start to couple the triangles and dimers with their 
neighbors, the spins are no longer free to move indepen- 
dently and they start to have a slower dynamics. In this 
case the magnetization obtained after a short time is dif- 
ferent from the equilibrium magnetization. 

A further confirmation of the very slow spin response 
to a magnetic field is given by the difference between the 
magnetization curve as a function of the time for the FC 
and the ZFC protocols, see Fig. 111 the following, we 
concentrate on the case T/J' = 10 -6 and h/J' = 10~ 2 . 
Let us first consider the ZFC magnetization. As long as 
h = 0, the magnetization vanishes and the chiralities are 
all opposite to the vectors e M and there are both triangles 
and dimers in the coupling pattern. When the external 
magnetic field is turned on, the magnetization reaches 
the value m sa 0.004 in about 20 Monte Carlo steps. This 
first very rapid regime corresponds to the spins moving 
toward the direction of the magnetic field. During this 
fast spin dynamics, the chirality does not play any impor- 
tant role in the dynamical process. Indeed, we have veri- 
fied that m ~ 0.004 is also the magnetization acquired in 
the first 20 Monte Carlo steps within a simulations where 
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FIG. 13: Long-time (upper panel) and short-time (lower 
panel) behavior of the magnetization per spin versus time 
at T/J' = 1CT 6 for a FC and a ZFC protocol. The external 
field has magnitude h/J' = 10 -2 . 



the chiralities were kept fixed. In the second regime, the 
magnetization increases very slowly up to its saturation 
value m ~ h/J'. This slow increase is due to the com- 
bined dynamics of spins and chiralities, in order to find 
a compromise to minimize the free energy. By contrast, 
within the FC protocol, the system is cooled down in the 
presence of the external magnetic field and therefore, at 
the final temperature, the magnetization has already a 
finite value m ~ 0.0085 and the chiralities arc no longer 
in the direction of the basis vectors. Then, by keeping 
the system at the desired temperature T < T m (h) , a very 
slow increase of the magnetization is observed, similar to 
the behavior of the second regime of ZFC. 

The pronounced difference between the ZFC and the 
FC protocols, and more generally an evident aging effect 
has already been emphasized in a related, but much sim- 
pler classical model containing two different Ising vari- 
ables^ Finding the same kind of behavior in a more 
realistic Hamiltonian demonstrates that this is a robust 
feature of models with competing variables and frustrat- 
ing interactions, even in the absence of an explicit disor- 
der like in ordinary spin-glass materials. 



V. CONCLUSIONS 

In this paper we have investigated the properties of 
the classical version of an effective Hamiltonian for the 
S = 1/2 Hcisenberg antiferromagnet on a trimerized 
kagome lattice using classical Monte Carlo methods. The 
analysis of this model shows many interesting features, 
possibly reminiscent of a spin-glass behavior at low tem- 
peratures. The fully frustrated character of the inter- 
action, due to the presence of two kinds of variables on 
each site, induces a freezing of the chiralities and, in the 
presence of an external magnetic field, the development 
of very large time scales for the spins at low temperature. 



Indeed, we have demonstrated that, by a partial order by 
disorder effect, for T < T* ps 0.05 J', the states with high- 
est entropy are selected by thermal fluctuations, and the 
chiralities freeze along particular directions. This pre- 
vents the system from developing any long-range corre- 
lations and completely splits the lattice into disconnected 
dimers and triangles, with antifcrromagnetically coupled 
spins. On each of these elementary entities the spins are 
free to rotate, as long as they add up to zero. 

When the spins are coupled to an external magnetic 
field and the temperature is sufficiently small, the freez- 
ing of the chiralities has a strong effect on the low- 
temperature spin properties. First of all, the equilibrium 
magnetization for small magnetic fields has a clear non- 
linear behavior, whereas for large enough h, we obtain 
m oc h. The actual magnetic field range for which the 
non-linear behavior is detected strongly depends upon 
the temperature. However, the non-linearity cannot be 
meaningfully casted into a divergent X3> m contrast to 
traditional spin glasses. Still, very strong non-linear 
effects in the magnetization curves are clearly present 
and point to a spin-glass reminiscent behavior. Sec- 
ondly, there is a clear difference between the FC and 
the ZFC measurements of the magnetization, typical of 
spin-glass models. A difference with ordinary spin-glass 
systems, representing an original aspect of our results is 
that the characteristic temperature T m below which the 
slow spin response takes place does not coincide with T* 
and strongly depends on the external magnetic field, i.e., 
T m = T m (h). On the other hand, the freezing transition 
of an ordinary spin glass corresponds to a thermodynamic 
transition and the freezing temperature is independent of 
the particular protocol. Nonetheless, the very slow spin 
response to an external magnetic field is typical of ordi- 
nary out-of-equilibrium systems, and here it is driven by 
the freezing of the chiralities. In this respect, the rich- 
ness of our model consists in having two distinct temper- 
atures, i.e., T* , which characterizes the entropic selection 
of chirality configurations, and T m (h), which character- 
izes the establishment of the slow spin dynamics. On the 
contrary, in the ordinary spin-glass models, there is only 
one characteristic temperature. 

So the results we have obtained on the classical version 
of the effective model of the S=l/2 Heisenberg model on 
the trimerized kagome clearly show that the interplay 
of two local degrees of freedom can lead to a very slow 
dynamics, and to a behavior which is in many respects 
reminiscent of the spin-glass phenomenology One point 
worth noticing at that stage is that this behavior is, as 
far as we can tell, independent of the dynamics chosen 
in the numerical simulations. In particular, we have per- 
formed a few simulations with simultaneous updates of 
spins and chiralities at the same site, and the results 
were unchanged. This is a clear improvement over the 
simplified model of Ref. I28L in which the entropic bar- 
riers were related to the dynamics. The present results 
suggest that the development of a very slow dynamics is 
actually a generic feature of models with two local de- 
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grees of freedom. 

Coming back to the experimental evidence in favor of 
a spin-glass behavior in SrCr8_ x Ga4+ x Oig, the present 
results arc very promising in several respects. First of 
all, the model was motivated by physical considerations 
of the exchange processes in this system, in which the 
kagome layers are trimcrized. Besides, the broken ergod- 
icity which is at the origin of the very slow dynamics is 
neither related to a specific dynamics nor to disorder but 
to the presence of two local degrees of freedom, a very 
natural consequence of the special topology of the trimcr- 
ized kagome lattice. In particular, although we started 
from spins 1/2, all half-integer spin models on the trimcr- 
ized kagome lattice lead to very similar effective models 
with two local degrees of freedom^ In addition, the char- 
acteristic temperature at which large time scales appear 
is much lower than the coupling constant, again in agree- 
ment with SrCrg_ x Ga4 +x Oig. Finally, the spins have a 
very anomalous response to magnetic field, with a clear 
difference between field-cooled and zero-field-cooled pro- 
tocols, while retaining a very fast dynamics at low tem- 
peratures. As far as we know, this is a unique feature of 
the present model, which could be related to the muon 
spin rotation observation that only a tiny part of the 
spins are actually frozen in SrCrs_ x Ga4 +x Oig below the 
spin-glass transitional 

There are clear differences however between our model 
and the properties of SrCr8_ x Ga4 +x Oig. An obvious one 
is the specific heat, which is quadratic below the spin- 
glass transition, whereas in our model it goes to a con- 
stant. This behavior is clearly a consequence of our ap- 
proximation to treat spin and chirality degrees of freedom 
as classical variables. Given the quasi-two dimensional- 
ity of the compound, one can speculate that a quantum 
treatment would lead to a T 2 behavior of the specific 
heat for the quadratic degrees of freedom, in agreement 
with experiments. This could actually be a natural ex- 
planation of this behavior, which is not the most com- 
mon behavior in spin glasses, which tend to have a lin- 
ear specific heat at low temperature. At that stage, the 
main discrepancy between the phenomenology described 
in the present paper and the experimental results ob- 
tained in SrCr8_ x Ga4+ x Oig is the observation of a clear 
phase transition in this compound with a divergent non- 
linear susceptibility. It is likely that a true phase tran- 
sition would require to include inter-layer coupling, like 
in disordered models of spin-glasses. Whether including 
such a coupling can produce a divergent non-linear sus- 
ceptibility is left for future investigation. We should also 
note that the effective model for spins 1/2 has specific 
features with respect to the effective models for larger 
half-integer spins. In particular, the identification of ex- 
act ground states does not seem to be possible for larger 
spins^ and the phenomenology of these models might 
turn out to be different in some respects, although we 
are confident that the presence of two degrees of freedom 
will still lead to an anomalous response of the spins to an 
external magnetic field and to the development of large 



time scales. 
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APPENDIX A: MONTE CARLO TECHNIQUE 

In this Appendix, we describe the Monte Carlo method 
used throughout our study. For simplicity, let us consider 
the Hamiltonian: 



(id) 



(Al) 



where &i are classical unit spin vectors on the sites of an 
TV-site lattice and Jy are magnetic couplings. 

In order to study the dynamics of this Hamiltonian, we 
have performed Monte Carlo simulations using the heat- 
bath algorithm to update the vector directions^ Each 
spin is assumed to be in contact with a heat bath and 
is immediately put into a local equilibrium with respect 
to the instantaneous effective field on it from the nearest 
neighbor spins: 



U)i 



(A2) 



where indicates the nearest neighbors of the site i 
and T is the temperature. 

The new direction of a given spin is determined by a 
polar angle 9 and an azimuthal angle <f> relative to Hi, 
through the following Boltzmann distribution: 



P(cos(9) = 



H 



Hcos 9 



2sinh_ff 



(A3) 



where H = \Hi\. The two angles 9 and <j> are chosen with 
the use of two random numbers Ri and i?2, uniformly 
distributed between and 1. Even though l|A3(l is inde- 
pendent of 0, it is useful to update this variable in order 
to increase the rate at which the phase space is sampled. 
We therefore chose a new </> with 

(j} = 2nR 2 . (A4) 

The polar angle 9 is instead obtained by 

cos0 = — ln[l + i2i(e 2i? - 1)]-1. (A5) 
H 

It is important to note that the heat bath method does 
not rely on an accept-reject approach. Instead, the new 
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configuration is already generated with the right proba- 
bility. 

When both spins and chiralities are present in the 
Hamiltonian, the same updating algorithm can be ap- 
plied to the chiralities. In our simulations, we consecu- 
tively chose a random site on the lattice, updated the spin 
and then the chirality on another random site. Quanti- 
ties like the energy or the magnetization where usually 
calculated after N such moves, i.e., after the system had 
undergone one Monte Carlo step per site on average. 

Because we are working with continuous degrees of 



freedom and not with discrete variables as in the Ising 
model for example, the heat bath algorithm is better 
suited than the conventional Metropolis method. In the 
Metropolis method, at low temperatures, one starts re- 
jecting many moves because the system prefers to make 
small fluctuations about local minima of the energy. 
Thus, it may be very difficult to reach equilibrium or 
even to know whether one has reached it. By contrast, 
the heat bath method provides an efficient way to reach 
the actual minimum of the energy because a new vector 
orientation is obtained every time, at any temperature. 
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